Non-equilibrium dynamics of coupled qubit-cavity arrays 
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We study the coherence and fluorescence properties of the coherently pumped and dissipative 
Jaynes-Cummings-Hubbard model describing polaritons in a coupled-cavity array. At weak hopping 
we find strong signatures of photon blockade similar to single-cavity systems. At strong hopping 
the state of the photons in the array depends on its size. While the photon blockade persists in a 
dimer consisting of two coupled cavities, a coherent state forms on an extended lattice, which can 
be described in terms of a semi-classical model. 



The realization of effective photon-photon interactions 
in various cavity QED systems has triggered an im- 
mense interest in using these light-matter systems for 
quantum computation and simulation. A central ele- 
ment of this work is the photon blockade [l| , where the 
presence of a single photon in a driven cavity prevents 
more photons from entering. Key experimental signa- 
tures of this effect are anti-resonant lineshapes in ho- 
modyne/heterodyne detection 0] as well as photon anti- 
bunching and the appearance of a dressed state MoUow 
triplet in resonance fluorescence spectra So far, these 
signatures have been observed in single cavity systems. 

With the unprecedented experimental control of recent 
experiments at hand, a key challenge in the study of 
coupled light-matter systems is the interplay of strong 
correlations and collective behavior in extended systems. 
Recent experimental progress includes the realization 
of a Tavis-Cummings nonlinearity with superconducting 
qubits in circuit QED 0, 0] , a BEC of weakly interact- 
ing exciton-polaritons in semiconductor microcavities Q 
and a non-equilibrium superradiant state using ultra-cold 
atoms in optical cavities Theoretical interest has fo- 
cused on a possible Mott insulator-superfluid transition 
of polaritons in coupled-cavity arrays as described by the 
Jaynes-Cummings-Hubbard model (JCHM), where each 
cavity is strongly coupled to a two-level system (2LS) and 
photons can hop between cavities [sl-fisl. However, with 
a few exceptions 3-13], most of these studies did not 



take into account the basic nature of quantum optical 
applications, drive and dissipation. 

In this paper, we study the coherently pumped JCHM 
including dissipation via spontaneous emission of the 2LS 
and cavity loss. We compare exact numerical simulations 
of a dimer of two coupled cavities with mean-field theory 
for an extended array. We calculate the pump frequency 
dependence of the photon field as measured in a homo- 
dyne/heterodyne detection scheme, the second-order co- 
herence function (photon statistics), and the fluorescence 
(emission) spectra. At weak hopping, we find strong sig- 
natures of photon blockade as observed in single-cavity 
systems. At strong hopping, the state of the photons 



in the driven-dissipative array depends on its size. In a 
dimer model consisting of two coupled cavities we find 
a photon blockade even at large hopping. For a lattice 
consisting of infinitely many coupled cavities the block- 
ade effect vanishes and a coherent photon state emerges, 
which can be described scmi-classically. We find that the 
crossover from weak to strong hopping is smooth as long 
as the array is pumped at the bottom of the lower polari- 
ton band. Higher pump frequencies may cause tunneling- 
induced bistabilities. 

The driven JCHM is described by the Hamiltonian 
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'0 



(1) 



Here, hj'^ denotes the Jaynes-Cummings model (JCM) 
at site i of the array. 



(2) 



Each 2LS is represented by a spin operator and coupled to 
the cavity photons with strength g. The second term in 
([!]) describes the hopping of photons to one of its z near- 
est neighbors at a rate J/ z (such that the total bandwidth 
is 2J). All cavities are subject to a coherent laser drive of 
strength / described by the third term in ([T]). Working 
in a frame that rotates at the pump (laser) frequency , 
we measure the energies of resonators (cavities) and 2LS 



(qubits) with respect to the pump, (5, 



(9) 



Dissipation is taken into account in the master equation 
for the density matrix 



p] + («;/2)i^N + (7/2)^K 



(3) 



with if [cr^ ] = Y.i (2(^i per. - tr+cr^ p - paja^ ) and a 
similar decay term for the photon operators. Here, k de- 
notes the cavity decay rate and 7 the spontaneous emis- 
sion rate of the 2LS. 

Recently, it was shown that the single driven and dis- 
sipative JCM (J = 0) for zero detuning (w^ = '-^q) yields 
surprisingly complex behavior 0, [20| . which we briefly 
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FIG. 1. Level scheme for the lowest excitations of the JCHM 
for (a) a single cavity, (b) a dimer model, (c) a coupled cav- 
ity array. Hopping between cavities gives (anti-)symmetric 
superpositions of photon states in (b) and photon bands 
e(k) = Jcos(k) in (c). The qubits Q are resonant with the 
lowest photon states: (a) the cavity mode C {ujq = ujr), (b) 



the symmetric superposition of photon states {uiq 



J), 



(c) the bottom of the photon band {ujq = ujr — J). This gives 
rise to dressed (polariton) states (LP, UP) in (a) and (b) and 
polariton bands in (c). In this paper the laser frequency uii 
is near the lowest excitations in the system (bottom of the 
polariton band). 



summarize here. The local Hamiltonian in ^ has eigen- 
states that are symmetric and anti-symmetric superposi- 
tions of excited 2LS and photon states, upper and lower 
polariton dressed states (Fig. [T^). Pumping weakly near 
the lower polariton, 8r ~ g, yields a resonance in the 
photon field 4> = \{a)\ with Lorentzian lineshape (dashed 
line in Fig. [2^). However, beyond linear response the ho- 
modyne signal turns into anti-resonant behavior and the 
Lorentzian develops a central dip of width ^ f / n (black 
line in Fig. [2^). This effect can be understood by re- 
stricting the Hilbert space to two states, the vacuum and 
the state with a single lower polariton (LP). The anti- 
resonance arises when this effective 2LS saturates. This 
effect is the semi-classical Ra bi sp litting corresponding to 
the dressing of dressed states 2l| which has recently been 
observed in circuit QED as a MoUow triplet in fluores- 
cence spectra . The 2LS approximation is appropriate 
as the nonlinearity J7off = £2 — 2ei = g{2 — a/2) of the 
JCM prevents higher states from being excited (ei,2 are 
the lowest energies in the Hilbert space sector with 1 or 
2 excitations, respectively). 

We now consider extended systems as described by 
the JCHM in Eq. ([J) with J ^ 0. The level scheme 
for the lowest excitations is shown in Fig. [U where we 
choose the qubits to be resonant with the symmetric pho- 
ton state of the dimer consisting of two coupled cavities 
(Fig. [T|d) or the bottom of the photon band in the in- 
finite array (Fig. [T};). At weak hopping we expect the 
coupled-cavity array to exhibit a similar blockade effect 
as a single cavity. At strong hopping the effective non- 
linearity in the spectrum to leading order in g ^ J is 
U^s = 2g{l- ^1 - l/{2Ns)^ for a ID chain. This only 

vanishes as the number of lattice sites Ns becomes large. 
Small arrays are thus expected to always show photon 
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FIG. 2. Photon field |0| and second-order coherence 
g'^-* for the dimer of two coupled cavities, (a) as a 
function of resonator detuning 5r/ff at hopping Jjg = 
(black), 0.02 (red), 0.05 (green) and drive strength fjg = 
0.005, and (b) as a function of hopping strength Jjg at drive 
//5 = 0.001 (black), 0.005 (red), 0.02 (blue) when the laser fre- 
quency is resonant with the lowest excitation (5^ = <7 + J , see 
Fig. [ij . The dashed line in (a) corresponds to a single cavity 
(J = 0) at very low drive //g = 0.001. The crosses mark 
points in (a) where 5r = (/-f J. Panel (c) shows g'^' as a func- 
tion of the hopping strength J/ g for the same drive strengths 
as in (b). All dissipation rates are k = 7 = 0.005(;. 



blockade, large ones only at small hopping. 

Figure [5] shows the photon field = and the 

second-order coherence 5'^-' = (a|a|aiai)/(a|ai)^ for a 
dimer model consisting of two coupled cavities as ob- 
tained from exact numerical evaluation of the master 
equation in ([3]). In Fig. [5^, the anti- resonance broad- 
ens when the hopping strength increases and shifts to 
larger values of detuning 5^/(7, i.e., smaller pump fre- 
quencies. The appearance of a second anti-resonance for 
J = 0.05g is associated with a two-polariton state with 
one polariton in each cavity. A similar effect has been re- 
ported for a three-site ring lattice in [l^ . The right-hand 
panel in Fig. [5] shows the crossover to the large hopping 
regime for several values of the pump strength. The anti- 
resonant feature remains for J ^ 5, i.e., the photon field 
decreases as a function of pump strength at fixed hop- 
ping. The second-order coherence g^^^ in Fig. [2] confirms 
this photon blockade picture. At weak drive (/ ^ C/off) 
photons in the dimer remain anti-bunched with g*-^) <C 1 
even at strong hopping. Note, that the dimer model is re- 
alizable with current state-of-the-art technology and has 
been shown to exhibit interesting new physics due to the 
coupling between cavities in various other contexts [22- 
23. 

We now increase the coordination number z and 
consider a large array of coupled cavities. Continuous 
bands form between the upper and lower polariton 
states (Fig. [T]:) and the full problem is then no longer 
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FIG. 3. Photon field |(j!>| (upper panels) and second-order 
coherence g'^' (lower panel) for the array as obtained from 
mean-field theory for the same parameters as in Fig. (2] 
Dashed lines in (a) and (b) correspond to the effective spin 
model in Eq. Q. The solid horizontal lines in (b) are the 
semiclassical asymptotes for Eq. ((5]). 



numerically tractable. For large z, the hopping term 
is approximated by ajaj « {(4)o.j + ^ {o\){o-j)- 

As all cavities are equivalent, we take the photon 
expectation values to be the same at each site of the 
array (f> = (ai). This is the same mean- field decoupling 
as in the corresponding equilibrium model Q. It was 
first used out-of-equilibrium to calculate the dynamics 
of the dissipative Bose- Hubbard model [l^. One may 
then find the steady state of this decoupled on-site 
problem numerically (Fig. [3|). As in the dimer model, 
the anti-resonance shifts and broadens with increasing 
hopping. The two sidepeaks of the anti-resonance be- 
come asymmetric and turn into a jump corresponding to 
a hopping-induced bistability. Similar bistable behavior 
in quantum optics is well known for a single strongly 
driven cavity [isl - fsH and for the Dicke model [H] . Here, 
a bistability develops at fixed pump strength when in- 
creasing the hopping rate J between cavities. Note, that 
these instabilities only appear when the array is pumped 
above the bottom of the lower polariton band. In Fig.[3)D 
we show the crossover to large hopping. We pump at 
the bottom of the polariton band, where no instabilities 
occur. The curves for different pump strengths cross 
and reorder, such that at strong hopping the photon 
field increases uniformly with drive strength. Thus, 
the anti-resonance feature vanishes. The second-order 
coherence function g*-^^ undergoes a smooth crossover 
from anti-bunching at weak hopping with g^^' w to a 
coherent photon state at strong hopping with g^^^ w 1 
even at vanishingly small drive strength. Thus, all 
signatures of the photon blockade effect are destroyed in 
a large array at strong hopping. 
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FIG. 4. Fluorescence spectra S{ijj) for the dimer (upper pan- 
els) and for the array (lower panels) at fixed drive strength (a) 
/ = O.OI5, (c) / = 0.0053 and fixed hopping (b,d) J = IO5. 
The laser frequency is resonant with the lowest polariton state 
{Sr = g + J, see Fig. [1}. Blue crosses and red circles in (c, d) 
correspond to the effective models at weak drive in Eq. 
and strong hopping in Eq. ([5]), respectively. The decay rates 
are k = 7 = O.OOlf? in (a, b) and k = 7 = O.OOSg in (c, d). 



Another important experimental signature is the fiu- 
orescence spectrum. We calculate the emission spec- 
trum as the Fourier transform S{uj) of the on-site auto- 
correlation function S{t) = {a^{t -\- T)a{t)) in the steady 
state of the system. Results are shown in Fig. |4l The 
spectrum of the dimer (Fig. |4]) exhibits a MoUow triplet 
of emission peaks at weak hopping at a; sa 0, ±-\/2/ aris- 
ing from transitions between the vacuum and the lower 
polariton state (LPi) dressed by the laser field. This 
agrees with recent circuit QED experiments in a single 
cavity and can be interpreted as a signature of photon 
blockade Q. For large hopping, the low-energy spectrum 
again exhibits a Mollow triplet, however with modified 
frequencies near w « ±2/, which is due to a shift of the 
weight of the photonic component in the lower polariton 
state LPi with increasing hopping. 

The lower panel in Fig. |4] shows the fiuorescence 
spectra for the array as obtained from mean-field theory. 
At small J, the lower polariton mode at a; « splits and 
the corresponding side-peaks move to higher frequencies 
as hopping increases. Despite the absence of the photon 
blockade the triplet structure survives at large hopping. 
This is because the resonant frequency of a single cavity 
is different from that of the array — it is detuned by 
J + g from the pump. In the limit of small drive /, where 
the photon field vanishes, the emission spectrum is hence 
due to off-resonance fluorescence of pumping a single 
2LS at g below its transition frequency, giving peaks at 
±g (Fig. HJi). Additional peaks with vanishing spectral 
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weight as J — >■ arise at higher energies uj ss ±2(7 
due to the emission from upper polariton states (for 
the dimer these modes are not visible on the scale shown) . 

Further insight into the behavior of the array at weak 
driving (/ <C C/cff, where UcS — >■ when J ^ g) can 
be gained through an effective model by generalizing the 
two-level dressed state approximation derived for a sin- 
gle cavity in Q to the JCHM. By projecting into this 
restricted Hilbert space, we obtain an effective Hamilto- 



z ^ 



.(4) 



which is valid when the number of excitations per 
cavity is small. The Pauli matrices {t+,t^,t^} de- 
scribe the transition between the vacuum and the 
lower polariton at energy t = 8r svy? 6 + Sq cos^ 9 
with e = tan-^{-2g/{5r - 5g))/2. Here, J = J sin^ 6 
and / = /sin(0) are effective hopping and pump 
strengths. The master equation for the effective system 
becomes p = —i[Heg,p] + Kl£\r^\ with the decay rate 
k = n%\T? Q -I- 7COS^6'. This model is equivalent to a 
dissipative Heisenberg XX m odel in a magnetic field 
Dissipative spin chains 
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with dissipation acting 
only on the boundary sites, have recently been con- 
sidered as examples of non-equilibrium quantum phase 
transitions. The current problem, however, differs from 
these in having dissipation on all sites. Dissipative spin 
chain problems also arise in the context of the Rydberg 



blockade [36|, |37[. Using the same mean-field approxi- 



mation as above produces an effective 2LS Hamiltonian 
-ffeff = (e/2)Tz -I- (/offT+ -I- h.c.) with an effective drive 
strength /cff = Jif,' + f and the self-consistent field 
if,' = (r^). The dashed lines in Fig. [3] correspond to 
the steady state of this model and agree well with the 
numerical treatment at small drive and hopping. Since 
UeS — )• for large J, the blockade picture breaks down 
with increasing hopping. The emission spectra resulting 
from this approximation are shown as blue crosses in 
Fig.H 



Pumping at the bottom of the polariton band (Fig. [T]:) 
we have 6r — g + J, such that the dependence on hopping 
disappears for the steady state in ([5]). Fig. [Sb shows 
that the corresponding asymptotes (solid horizontal 
lines) are indeed approached by the mean-field numerics 
for J — > oo. The low-energy spectra S{ijj) in this 
approximation are found by linearizing equations ([5]) 
about their steady state. Fluctuations on different sites 
are uncorrelated in ([5]) and the 1/z term can thus be 
neglected. This introduces a J-dependence into the 
spectrum. Fig. U shows that the predictions of this model 
match the full numerics of the lattice mean-field theory 
everywhere except at very small hopping (deviations not 
visible on the scale shown). The equations of motion 
in ([5]) are the same as those obtained from a Dicke 
model where the single site expectation values in ([S]) 
are replaced by collective spin variables. However, the 
spectra differ from the usual Dicke model in describing 
on-site rather than collective fluctuations. 



In summary, we have investigated the crossover from 
weak to strong hopping in the coherently driven dissipa- 
tive JCHM. In the weak hopping limit, blockade physics 
survives and both the extended and two-site arrays be- 
have similarly. At large hopping, blockade effects survive 
in the two-site system, while the infinite system becomes 
scmiclassical and can show bistability for certain pump 
frequencies. Our work thus provides a strong motiva- 
tion for further experimental and theoretical studies of 
coupled cavity arrays. 
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The strong hopping regime of the full lattice can be un- 
derstood by relating the JCHM to a pumped dissipative 
semi-classical model. Rather than restricting the Hilbert 
space to low excitation numbers, we let the photons oc- 
cupy a coherent state and factorize products of spin and 
photon expectation values. The on-site equations of mo- 
tion are 



0, = -liSrcP, + / + gsr - (J/z) ^ 0,) - («/2)(/., 



■h = -iSqS, + 2ig(l),sl - (7/2)s.; 

st = -^g{<p^s-* - (j)*s^) - 7 (sf + 1/2) , 



(5) 



where st = (ct 



111% 



(erf) and 



{ai). In the steady 



state we assume that all sites are identical 
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